Scale-invariant Truncated Levy Process 
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We develop a scale-invariant truncated Levy (STL) process to describe physical systems character- 
ized by correlated stochastic variables. The STL process exhibits Levy stability for the probability 
density, and hence shows scaling properties (as observed in empirical data); it has the advantage 
that all moments are finite (and so accounts for the empirical scaling of the moments). To test the 
potential utility of the STL process, we analyze financial data. 



In recent years, the Levy process [l| has been proposed 
to describe the statistical properties of a variety of com- 
plex phenomena (2|-g| . The Levy process is characterized 
by "fat tails" (power law), and display scaling behavior 
similar to that observed in a wide range of empirical data. 
However, the application of the Levy process to empiri- 
cal data is limited because it is characterized by infinite 
second and higher moments, while empirical data have 
finite moments. 

Truncated Levy (TL) processes are defined to have a 
Levy probability density function (PDF) in the central 
regime, truncated by a function decaying faster than a 
Levy distribution in the tails [Q. The TL process is 
introduced to account for the finite moments observed 
for empirical data ||ll[. However, the TL process (with 
either abrupt Effl or smooth [Q truncation) has limita- 
tions when applied to empirical data, (i) The TL process 
is introduced for independent and identically distributed 
(i.i.d.) stochastic variables, while variables describing 
many physical systems are not i.i.d. — e.g. there are 
correlations, (ii) The PDF of the TL process tends to 
the Gaussian distribution (according to the central limit 
theorem), and hence does not exhibit scale invariance; 
PDFs for a variety of complex systems, however, are of- 
ten characterized by regions of scale-invariant behavior. 
(iii) The time scale above which the Levy profile becomes 
Gaussian depends on the size of the truncation cutoff (or 
the standard deviation) |10 IJ]; to mimic the Levy type 
scale invariant behavior observed for data, the TL pro- 
cess must be defined with a standard deviation larger 
than the one observed for the data [see caption to Fig.l]. 

Here we introduce a type of stochastic process which 
we call the scale-invariant truncated Levy (STL) pro- 
cess. Stochastic variables z in the STL are gener- 
ated by the symmetrical probability function f{z) = 
Ae~^'^l'^|z]_^~", where < a < 2 ||l|]. The exponential 
prefactor |12] ensures a smooth truncation of the Levy 
distribution, where the parameter (3 can take any posi- 
tive value [O, A~^ is related to the size of the truncation 
cutoff, and A is a measure of the "spread" in the central 
region. 

From the probability function /(z), we calculate 
the characteristic function 4>{k) = exp[— J_ dz(l — 
g-ifcz-jj^-^)] ^ rpj^g pQp 75^^(2) is the Fourier trans- 



form of 0(fc) m 



V{z) ^ — I 4>{k)e^^^dk, 



(1) 



Since f{z) w ^|z|^-'-^" for small values of z, V{z) has a 
Levy profile in the central part. To maintain scale invari- 
ance for V{z) in the entire range including the tails p6[ , 
we introduce the scaling transformations 



^At = (Ai)Mi, AAt = (At)-^''/"Ai 



(2) 



where Ai is the time scale and e can take any positive 
value. Under these transformations, the PDF V{z) = 
VAt{z) scales as the Levy stable distribution: 
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With the transformations of Eqs. (||) and (||), we obtain 
a process with controlled dynamical properties — 'pAt(-z) 
for any value of At can be calculated from the PDF at 
any chosen At (e.g. At = 1) Q. 

Although the PDF V\t{z) exhibits scale invariant 
properties identical to the Levy stable distribution, the 
process defined by Eqs. (|^) and (||) is different. While 
the Levy process is defined for i.i.d. variables the STL 
process is characterized by correlated stochastic variables 
— the STL is a non-i.i.d. type process. To demonstrate 
this, we consider the scaling of the second moment u"^ , 
determined as the second derivative of 4>{k) at small val- 
ues of /c H: 
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where a\ is the initial standard deviation for At = 1. 
The second equality on the right hand side follows from 
the transformations of Eq. (0) . For an appropriate choice 
of e/a {^ 0.5), the scaling relation (Q) indicates the pres- 
ence of correlations that can be positive (or negative). 
In addition, the STL process exhibits scaling not only 
for the second moment but for all higher moments: 

< Izl" >= I dz |z|" VAt{z) = Ai"'/" < |zi|" > . (5) 



Hence, the STL is a process for which the PDF VAt{z), 
the second moment a^, and aU higher moments < \z\" > 
scale with the same scahng exponent e/a. 

Often with empirical data, we observe several differ- 
ent scaling regimes. To account for a crossover at given 
time scale ( Ai) x , we introduce different scale invariant 
transformations from the type of Eq. (0) for two different 
regimes of time scales: 



^^* " \ (At)-^^^/"Ax 



AAt = 



{AtY'A, 



1 < At < (At)> 
Ai> (At)>, 



1 < Ai < (Af)> 
At > ( Ai) > 



(6a) 



(6b) 



Here a, Ai and Ai are free parameters, chosen to fit 
T^Atiz) at the initial time scale At = 1. Continuity 
of the PDF and the moments between the two scaling 
regimes is ensured by continuity in the values of A and 
A: from Eqs. (|^ and (|^) we find Ay, = {Aty^-'^Ai 



and A^ 



(At) 



I3{ei-e2)/a 



Ai. 



To exemplify the features of the STL process for a 
broad range of time scales, we need sufficiently large data 
sets. Such a large data set is the S'&P500 stock index 
over the 12 year period Jan '84-Dec '95. The price fluctu- 
ations z of this index are the stochastic variable analyzed. 
In particular, we focus on the scaling behavior of several 
statistical characteristics: (1) the second and higher mo- 
ments, (2) the probability of return to the origin 7'At(0), 
and (3) the PDF PAt{z). For simplicity we set /3 = 1. 

We make three empirical observations, (i) Experimen- 
tal results for the standard deviation as a function of At 
show two different scaling regimes with a crossover at 
(At)x ~ 30 min ||ll| [Fig.l]. The regime at small time 
scales is characterized by slope 0.7, indicating the pres- 
ence of positive correlations in the price fluctuations z 
( "superdiffusive" regime). The second regime has slope 
0.5, indicating absence of correlations ("normal diffu- 
sion" regime). Therefore the fluctuations in the S'&P500 
index cannot be described by an i.i.d. stochastic process, 
such as the Levy or the TL process, (ii) The probability 
of return to the origin ■pAt(O), however, exhibits a Levy 
type of scaling for more than three decades [Fig. 2]. Such 
scaling for 7-'At(0) therefore indicates Levy scale invari- 
ance of the central part of the probability density, (iii) 
The scaling exponent of 7-'At(0) is identical to the expo- 
nent of the standard deviation in the first scaling regime. 
However, the crossover in the scaling of the standard de- 
viation is not followed by a change in the slope of VAtiO)- 

To account for the first empirical observations, we in- 
troduce a stochastic process with two different regimes: 
(a) a STL regime with Aai = {AtyAi and Aa* = 
(At)~'^/"Ai, to account for the superdiffusive behavior 
a ex (At)'=/" (Eq. ||) at short time scales At < (At)x 
[Fig.l]; and (b) a regime with breakdown of scaling de- 
fined by Aa* = Ax = const and Aai = (At)Ax for 
At > (At)x to account for the normal diffusive behavior 



a ex (At)i/2 (Eq^ (1) aj^fj Fig.l). This breakdown allows 
for a transition from a non-i.i.d. STL process to an i.i.d. 
TL process. 

The STL process in the regime At < (At)x accounts 
for the second empirical observation, the identical scal- 
ing exponent (e/a) experimentally observed for both the 
standard deviation a (Eq. 0) and the probability of return 
to the origin 7-'At(0) (Eq. ]^ and Fig. 2). From fitting the 
initial probability distribution 'Pi(z), we obtain a — 1.43. 
Since empirically the standard deviation scales with ex- 
ponent e/a — 0.7, we find that e — 1 for this process. 

Third, we find that the theoretical prediction for the 
STL process with a scaling breakdown is in good agree- 
ment with the empirical result for 7-'At(0) for more than 
three decades [Fig. 2]. We note that the transition at 
(At)x ~ 30 from STL (non-i.i.d.) process to a TL (i.i.d.) 
process in the scaling of a [Fig.l], does not imply a 
sharp transition in the scaling of 'PAi(O) from a Levy 
to Gaussian behavior [Fig. 2]. The reason is that the STL 
scaling regime (Eq. (|)), 'PAt(O) exhibits Levy scahng 
behavior (Eq. (||)) up to (Af)x ~ 30. In this scaling 
regime, a increases superdiffusively with exponent 0.7, 
that is much faster than 0.5 for an i.i.d. process. At the 
crossover scale (At)x , the standard deviation reaches the 
value CTx = (At)OVi. The value of CTx = (At)OVTL 
can be also related to an i.i.d. TL process with initial 
standard deviation axL > ci [Fig.l]. According to the 
central limit theorem, an i.i.d. TL process asymptoti- 
cally converges to a Gaussian process. Thus while in the 
short time regime (small At) the price fluctuation z over 
time At is a sum of correlated stochastic variables, in 
the asymptotic regime (large At), z can be treated as 
a sum of newly-defined independent stochastic variables 
with standard deviation cttl ■ Since such a Gaussian pro- 
cess is defined with large initial standard deviation axL, 
the transition from the Levy to the Gaussian behavior 
is delayed [Fig. 2]. The time scale (At)^ of this transi- 
tion can be calculated by equating the return probability 
'Pa* (0) for the Levy and Gaussian distributions Q . We 
find that (At)^ = S(At)x, where S « 70 [Fig.2]. Such 
a relation is interesting, since it explicitly connects the 
crossover from the Levy to Gaussian with the crossover 
from non-i.i.d. to i.i.d. process. 

Finally, we compare the empirical distributions of the 
price increments z of the S'&P500 index for different time 
scales At with the shape of the distributions obtained an- 
alytically [Fig. 3]. Good agreement between data and the 
theoretical distributions is observed both for the central 
part and for the tails. At small time scales, the scale- 
invariant behavior of VAtiz) is maintained in the entire 
range (Levy for the central profile, and exponential in the 
tails) due to the scaling transformations of the STL pro- 
cess (Eq. H). The crossover to an i.i.d. TL process at large 
time scales ensures a smooth transition to a Gaussian-like 
profile. We find that the proposed mechanism of a STL 
process, with breakdown, provides a reliable control of 
the dynamical properties of the PDF. 

In our analysis, we have considered the price fluctua- 



tions of the S'&PSOO index as the stochastic variable z. 
The choice of stochastic variable depends on the type of 
the stochastic process: e.g., for an additive process one 
considers increments, while for multiplicative processes 
the appropriate choice is relative increments. In finance, 
it is traditionally assumed that economic indicators arise 
from a multiplicative process, and correspondingly the 
preferred quantity to analyze is the rate of return or the 
difference in the natural logarithm of price. The addi- 
tive and multiplicative processes are related for high fre- 
quency data (small At) and short period of analysis, so 
the use of price fluctuations or rates of return lead to 
similar results. We find that even for low frequency data 
(large Ai) and for long period of analysis (up to 12 years), 
the results for the PDF and the standard deviation re- 
main similar for both the price fluctuations and the rates 
of return [Fig. 4]. 

We have proposed a stochastic process that even in the 
presence of correlations among the stochastic variables 
exhibits a Levy stability for the PDF. The STL process 
is characterized by identical scaling exponents for both 
the moments and the PDF. The STL process provides an 
unified dynamical picture to describe different statistical 
properties, and can be generalized for situations when 
the moments and the PDF exhibit different scaling be- 
havior. The STL process can be utilized — as we show 
in the case for financial data — not only for processes 
with a single scaling regime but also for physical systems 
with different regimes of scaling behavior. 
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with changing values of /3: (i) for /3 < 1 we have stretched 
exponential tails; (ii) for fl = 1 we have an exponential 
truncation (the case of WM)\ (iii) /3 > 1 we have pro- 
cess faster than exponential; and (iv) for 13 = 2 we have 
Gaussian behavior. 

B. V. Gnedenko and A. N. Kolmogorov, Limit Dis- 
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(Addison- Wesley, Cambridge, MA, 1954). 
In the case of TL process, Levy-type scaling transforma- 
tions for the stochastic variable {z} and the parameter 
A — z = {At)~^^°'zi and A = AtAi — ensure scale 
invariant behavior only for the central profile of 'P(z). 
Note that the STL process characterized by given a can 
scale with any scaling exponent e/a in contrast to the 
Levy stable process which scales with the scaling expo- 
nent 1/a. The parameter e controls the dynamics of the 
process — probability distributions characterized with 
the same a can exhibit different scaling behavior for dif- 
ferent values of e. E.g. for e = 1 and A = under the 
transformations of Eqs. (0) the probability density Viz) 
scales as the Levy stable process. 

We obtain the following analytic expression: B — 
[V2^c^lrAt=l(0)]^"/(2-"^ where C is the Levy PDF [Q. 
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FIG. 1. The S'&P500 index shows two regimes of scal- 
ing behavior for the standard deviation a. The correlated 
(superdiffusive) regime at small At corresponds to the STL 
process with slope e/a = 0.7. To account for the crossover to 
uncorrelated (normal diffusion) regime, we introduce a break- 
down in the scaling for the STL process: Aa* = Ax = const 
and AAt = (At) Ay for At > (At)x. The breakdown in the 
STL is equivalent to a transition to a TL process at large time 
scales. This TL process corresponds to an initial o^tl larger 
than the empirical ai. This is the reason for the delay (at 
time scale {At)s ~ 10'^) in the transition from Levy to Gaus- 
sian behavior observed for 'PAt(O) [Fig. 2]. Note, that the TL 
process with an initial standard deviation cri (as observed in 
the data) would exhibit for 'PAt(O) a transition from Levy to 
Gaussian at shorter time scales. 
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FIG. 2. S&iP500 data for the probability of return to the 
origin pAt(O) mimics Levy scaling for more than 3 decades 
in At. The slope and the intercept of the straight line which 
represents the scaling of the Levy distribution are determined 
from the parameters a — 1.43 and Ai = 0.0014 by fitting the 
initial PDF Vi{z) for the S&PSOO index. From the same fit, 
we obtain Ai — 0.7. These initial parameters are used to 
define the STL process. As expected, the STL process fol- 
lows the Levy scaling for ■PAt(O) at all time scales. The TL 
process (with cti = 0.07, identical to the empirical value) ex- 
hibits a transition at short time scales to the Gaussian process 
(with the same value of cri), in disagreement with the data. 
The STL process with a breakdown at (At)x, however, is in 
agreement with the data and explains the delayed transition 
(at (Ai)s ~ 10"^) to the Gaussian observed in the data. 



10" 



N 




10" 



10" 



FIG. 3. SkP500 probability distributions PAt(z) of price 
fluctuations z for different time scales Ai. With solid lines, 
we show the PDF of the STL, with breakdown process for 
the same time scales and parameters used in Fig. 2. Good 
agreement between data and the theoretical PDFs is observed 
for the central part. To reproduce better the experimentally 
observed change in slope of the far tails we use a = 1.43, 
Ai — 0.0028, and Ai — 2.6. The shape of VAtiz) changes as 
a function of Ai from exponential-like (for small Ai — STL 
non-i.i.d. regime) to Gaussian-like profile of the tails (for large 
Ai — TL i.i.d. regime). 
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FIG. 4. For the S'&PSOO index (Jan '84- Dec '95), we show 
how the standard deviation of price fiuctuations and rates of 
return completely overlap, after rescaling by the average price 
for that period. Note, that the standard deviation is calcu- 
lated for a long period of time (12 years) and for data at either 
high frequency (1 min) or low frequency (1 trading week). In 
(b) and (c) we show how the probability distributions of the 
same stochastic variables collapse after rescaling by the same 
average price. The collapse is obtained again for both low and 
high frequency data. We also subdivide the 12 year period in 
two equal subintervals, and find the same collapse. 



